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Abstract 



We compare the numerical predictions for heavy quark production in high 
energy hadron cohisions of the conventional QCD parton model and the kx- 
factorization approach (semihard theory). The total production cross sections, 
one-particle rapidity and pT distributions as well as several two-particle correla- 
tions are considered. Some of them can help to estimate QCD scale. The difference 
in the predictions of the two approaches is not very large, while the shapes of the 
distributions are slightly different. 
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1 Introduction 



The investigation of heavy quark production in high energy hadron colhsions is an 
important method for studying the quark-gluon structure of hadrons. Reahstic estimates 
of the cross section of heavy quark production as well as their correlations are necessary 
in order to plan experiments on existing and future accelerators or in cosmic ray physics. 

The description of hard interactions in hadron collisions within the framework of 
QCD is possible only with the help of some phenomenology, which reduces the hadron- 
hadron interaction to the parton-parton one via the formalism of the hadron structure 
functions. The cross sections of hard processes in hadron-hadron interactions can be 
written as the convolutions squared matrix elements of the sub-process calculated within 
the framework of QCD, with the parton distributions in the colliding hadrons. 

The most popular and technically simplest approach is the so-called QCD coUinear 
approximation, or parton model (PM). In this model all particles involved are assumed to 
be on mass shell, carrying only longitudinal momenta, and the cross section is averaged 
over two transverse polarizations of the incident gluons. The virtualities of the initial 
partons are taken into account only through their structure functions. The cross sections 
of QCD subprocess are calculated usually in the leading order (LO), as well as in the 
next to leading order (NLO) |^, 0, ^. The transverse momenta of the incident 
partons are neglected in the QCD matrix elements. This is the direct analogy of the 
Weizsaecker- Williams approximation in QED. It allows to describe reasonably the 
experimental data on the total cross sections and one-particle distributions of produced 
heavy flavour^, however it can not reproduce, say, the azimuthal correlations ^ of two 
heavy quarks and the distributions of total transverse momentum of heavy quarks pairs 
1^, which are determined by the transverse momenta of the incident partons. 

There is an attempt to incorporate the transverse momenta of the incident partons 
by a random shift of these momenta {kx kick) according to certain exponential 
distributions. This allows to describe quantitatively the two-particle correlations 0, 
but it creates the problems in the simultaneous description of one-particle longitudinal 
and transverse momentum distributions j^. At the same time this procedure has no 
serious theoretical background. While the shift of the order of Aqcd ((^t) ~ 300 
MeV) looks to be reasonable as having an possible origin in confinement forces at large 



distances, the values (kx) ~ 1 GeV 0, or even 3-4 GeV ||T0|, ^ should be explained in 
terms of the perturbative QCD. 

Another possibility to incorporate the incident parton transverse momenta is re- 
ferred to as fcj'-factorization approach [|T^, |T^, n, |T^, [T^, |T^, or the theory of semihard 



interactions [|1^, [T^, ^ ^ |2^, ^ . Here the Feynman diagrams are calculated tak- 
ing account of the virtualities and of all possible polarizations of the incident partons. 
In the small x domain there are no grounds to neglect the transverse momenta of the 
gluons, qiT and q2T, in comparison with the quark mass and transverse momenta, Pix- 
Moreover, at very high energies and very high the main contribution to the cross 



^New FNAL data Q on beauty production are 2-3 times higher than the NLO parton model 
predictions. 
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sections comes from the region of qiT ~ Pit or g2T ~ Pit, see Sect. 4 for details. The 
QCD matrix elements of the sub-processes are rather complicated in such an approach. 
We have calculated them in the LO. On the other hand, the multiple emission of soft 
gluons is included here. That is why the question arises as to which approach is more 
constructive. 

The predictions of all phenomenological approaches rely on the quark and gluon 
structure functions. The last ones are more or less known experimentally from the data 
of HERA, but unknown at very small values of Bjorken variable x < lO""^. However it is 
just the region that dominates in the heavy quark production at high energies^. A more 
serious problem is probably related to the fact that the NLO parton model formalism 
is based on the conventional structure functions whereas the fc^-factorization approach 
uses so-called unintegrated distributions which, at the moment, are known with the 
accuracy not good enough. 

In Sect. 2 we shortly discuss the conventional NLO parton model and show some 
numerical results together with the experimental data. The main formalism of the 
fcr-factorization approach is presented in Sect. 3. 



In Sect. 4 we present a comparison of results obtained with the help of kx- 
factorization and the parton model. The main goal of this Section is to demonstrate the 
differences in the qualitative predictions coming from the matrix elements. To simplify 
the calculations and to avoid various additional dependences we have used a "toy" gluon 
distribution which has only a reasonable qualitative behaviour and a fixed value of as- 



More accurate comparison ||29| between the conventional parton model results and 
the fc-r-factorization approach is given in Sect. 5 for the experimentally measured quanti- 
ties. Here we have used more realistic gluon distribution GRV94 |^0| which is compatible 
with the most recent data, see discussion in Ref. Predictions of the fc^^-factorization 

approach for several heavy quark correlations (some of them have never been discussed 
before) are presented in Sect. 6. 



2 Conventional parton model approach 

The conventional parton model expression for the calculation of heavy quark hadropro- 
duction cross sections has the factorized form 

a{ab QQ) = Y1 j dxidxjGa/iixi, iJF)Gb/jixj, fiF)d-{ij QQ) , (1) 

where Ga/i{xi, ftp) and Gb/j{xj, fip) are the structure functions of partons i and j in 
the colliding hadrons a and b, fip is the factorization scale (i.e. the value of the order 
of the maximal virtuality of incident partons) and cr{ij QQ) is the cross section of 

^For example, in the case of charm production, nic — 1.4GeV, at LHC, y/s — 14 TeV, the product 
X1X2 of two gluons (both xi and X2 are the integral variables) is equal to 4 • 10~* and a ppl icapability 
of the existing structure functions seems not to be clear at so small x, see discussion in ||25| . Another 
problem at very small x can be connected to gluon shadowing jl^, 
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the subprocess which is calculated in perturbative QCD. The latter cross section can be 
written as a sum of LO and NLO contributions, 

Q 

where fiR is the renormalization scale and fj;f as well as f^j^ and f^j'' depend only on 
the single variable 

ArriQ 

P ~ Z ; "S = XiXjSab ■ (3) 

(In the factor ln(yU^/mQ) we assume fiR = fip following In the case of different 
values of fin and fip, which is preferable for the description of the experimental data 
0, Eq. (2) becomes more complicated.) 

The expression (1) corresponds to the process shown schematically in Fig. 1 with 
QiT = QjT = 0. The main contribution to the cross section at small x is known to come 
from gluons, i = j = g. 




Fig. 1. Heavy quark production in hadron-hadron collision. The LO parton model 
corresponds to the case when qu = g2i = 0. 



The principal uncertainties of any numerical QCD calculation are the consequences 
of unknown scales Q, /iF and and the heavy quark mass, mq. Both scales (some- 
times they are assumed to be equal) are to be of the order of hardness of the treated 

^These uncertainties have to disappear after all high order contributions are summed up. There is 
the opinion that the strong (weak) scale dependence in LO or NLO means the large (small) contribution 
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process, however which value is better to take, mg, rriT = y^^q + Pt ■5, remains 
to be unknown. In principle, the uncertainties should not be essential because of the 
logarithmical dependence on these parameters. Unfortunately the masses of c and b 
quarks are not large enough and it leads to numerically large uncertainties (see e.g. 0) 
at the existing energies for fixed target experiments. 

Usually in the parton model the values 

fiF = f^R = mq (4) 
are used for the total cross sections and 



fiF = fiR = rriT = ^ml+pl (5) 

for the one-particle distributions 0. 

Here we calculate the total cross sections of heavy quark production as the integrals 
over their px distrubutions, i.e. with the scales (5). 

Both in the parton model and in the /cT^-factorization approach we take 

rric = 1.4 GeV, rub = 4.6 GeV , (6) 



for the values of short- distance perturbative quark masses p3| , |34 . 

Another principal problem of the parton model is the coUinear approximation. The 
transverse momenta of the incident partons, QiT and Qjt are assumed to be zero, and 
their virtualities are accounted for through the structure functions only; the cross sec- 
tion a{ij — > QQ) is assumed to be independent of these virtualities. Naturally, this 
approximation significantly simplifies the calculations. 

The NLO parton model calculations of the total cross sections of cc and bb produc- 
tion, as functions of the beam energy, for it~N and p — N collisions can be found in 
0. These results depend strongly (by the factor of several times) on the numerical 
values of quark masses as well as on the both scales, fip and fiR. Although there is a 
contradiction in some experimental data generally (at fixed target energy) they are in 
agreement with NLO parton model predictions. 

The NLO corrections to one-particle distributions lead only to renormalization of 
the LO results, practically without changing in the shapes of the distributions |3^, |[. 
It means that instead of the more complicated NLO calculation of the px or rapidity 
distributions it is enough to calculate them in LO, and to multiply then by K-factor 

which can be taken, say, from the results on the total production cross sections. The 
comparison of LO + NLO calculations with the LO multiplied by the K-factor is pre- 



sented in Fig. 2 taken from Ref. |35 



of high order diagrams, but it is not the case. The strong or weak scale dependence in LO or NLO 
indicates only the same dependence for the higher orders, which could be numerically small or large at 
some particular fixed scale. 
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Fig. 2. The calculated pj-distributions f or pp ^ b + X at ^/s = 1.8 TeV, for different 
values of rapidity for LO+NLO QCD parton model and for LO contribution multiplied 

by K-f actor, Eq. (7) equal to 2.5. 



The values of K-factors and their energy and scale dependences for the various sets 
of structure functions were calculated in Refs. |^ These dependences are similar 
for the modern sets of structure functions. 



The experimental data of WA92 PBI and E769 PP[ for Xi?-distributions of D-mesons 



produced in vrA^ interactions are in agreement with the parton model distributions for 
bare quarks, as one can see in Fig. 3 taken from [Q. It means that the fragmenta- 
tion processes are not important here, or they are compensated by, say, recombination 
processes. The shape of XiT'-distributions does not depend practically on the mass of 
c-quark. 

It was claimed in that the data on one-particle py-distributions, both obtained by 
fixed target experiments and by hadronic colliders (DO, CDF and UAl collaborations) 
for the case of beauty production can be described by the NLO parton model using the 
variation of both scales, fiR and fip- However new FNAL data on beauty production 
are 2-3 times higher than the theoretical curve 

More serious problem appears when we consider the correlations of produced heavy 
quarks. Let as start from so-called azimuthal correlation. The azimuthal angle is 
defined as an opening angle between two produced heavy quarks, projected onto a 
plane perpendicular to the beam and defined as the xy-plane. In the LO parton model 
the sum of the produced heavy quark momenta projected onto this plane is exactly zero. 
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Fig. 3. Experimental xp distributions for D mesons, compared to the NLO parton 

model prediction for charm quarks. 



and the angle between them is always 180°. In the case of the NLO parton model the (j) 
distribution is non-trivial [§] , however the predicted distribution (without including the 
kT kick) is narrower in comparison with the fixed target data 0. 

The investigation of such distributions is very important. The relative LO and NLO 
contributions of the parton model depend on the unknown scales, hf and hr. In the 
one-particle distributions, the sum of these contributions practically coinsides with the 
LO contribution multiplied by the /T-factor, as it is shown in Fig. 2. Therefore we 
are unable to control separately the magnitudes of the LO and NLO contributions. 
Contrary, in the case of azimuthal correlations all difference from the simple 5(0 — vr) 
distribution comes from the NLO contribution. 

The experimental data on azimuthal correlations are claimed (see [|D, ^ and Refs. 
therein) to be in disagreement with the NLO predictions, for the cases of charm pair 
hadro- and photoproduction at fixed target energies. The level of the disagreement can 
be seen in Fig. 4 (solid histograms) taken from Ref. 0. These data can be described 0, 
assuming the comparatively large intrinsic transverse momenta of incoming partons [kp 
kick). For each event, in the longitudinal centre-of-mass frame of the heavy quark pair, 
the QQ system is boosted to rest. Then a second transverse boost is performed, which 
gives the pair a transverse momentum equal to pxiQQ) = ^t(I) + ^r(2); kx^l) and 
A;r(2) are the transverse momenta of the incoming partons, which are chosen randomly, 
with their moduli distributed according to 

-^ = ^exp(-4/(4)). 



N dkl {kl) 

The dashed and dotted histograms in Fig. 4 correspond to the NLO parton model 
prediction, supplemented with the kx kick with = 0.5 GeV^ and {k?j.) = 1 GeV^, 
respectively. We see that with {k^) = 1 GeV^ it is possible to describe the data. 

However the large intrinsic transverse momentum significantly changes one-particle 
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4. Azimuthal correlation for charm production in vrA^ collisions: NLO parton 
model and kick calculations versus the WA75 and WA92 data. 



Pr-distributions of heavy flavour hadrons, which were earlier in good agreement with 
the data. The solid curves in Fig. 5 taken from |^ represent the NLO parton model 
predictions for charm quarks pT-distributions which are in agreement with the data. 
The effect of the kick results in a hardening of the spectrum. On the other hand, 
by combining the kr kick with (k^) = 1 GeV^ and the Peterson fragmentation |[42|| , 
the theoretical predictions slightly undershoot the data (dot-dashed curves). For the 
predictions at higher energy the values (k^) = 3-4 GeV^ [jlO|, |ll|] were used without any 
connection with perturbative QCD diagrams. 

The kx kick affects the xjT'-distributions of the produced c-quarks very weakly, Fig. 3, 
and after accounting the fragmentation these distributions become too soft. 

The general phenomenological expression for the heavy quark production can be 
written^ as a convolution of the initial transverse momenta distributions, I{qiT) and 
I{<12t)i with squared modulo of the matrix element, 

CTQCDiQQ) OC j d?qiTd?q2Tl{<llT)I{<l2T)\M{qiT,q2T,PlT,P2T)\^ ■ (9) 

Here there are two possibilities: 

i) the typical gluon transverse momenta are much smaller than the transverse momenta 
of produced heavy quarks, qiT -C PiT, or 

ii) all transverse momenta are of the same order, qi^ ~ piT- 

In the first case one can replace both initial distributions I{qiT) by 5-functions. It 
reduces the expression to the simple form (coUinear approximation): 

(JcoiiiQQ) OC |M(0,0,Pit,P2t)|' (10) 

in total agreement with Weizsaecker- Williams approximation. 
^We omit here for the simphcity ah unimportant factors. 
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Fig. 5. Charm p\. distribution measured by WA92 and E769, compared to the NLO 
parton model predictions, with and without the non-perturbative effects. 



In the second case we can not a priory expect good results from the Weizsaecker- 
Williams approximation. However it gives the reasonable numerical results in some 
cases. 

The kT kick |]^ discussed above effectively accounts for the transverse momenta of 
incident partons. It is based on the expression which symboUically reads 

(TkickiQQ) oc /(gir)/(g2T)® |M(0,0,piT,P2T)r (11) 

that has no good theoretical background. 

The main difference from the general QCD expression, Eq. (^, is that due to the 
absence of qiT in the matrix element the value {k"^) has to be differently taken for 
different processes (say, for heavy flavour production with comparatively small or large 
Pt)- The reason is that the functions I {qiT) in general QCD expression decrease for large 
Qj^j^ as a weak power (see next Sect.), i.e. comparatively slowly, and the q^rp dependence 
of the matrix element is more important. 

In the last one the corrections of the order of qfrp/ yu^, where /i^ is the QCD scale, are 
small enough when << 1 and they start to suppress the matrix element value 

when glr/^i^ ~ 1. Thus the essential values of the q^j, depend on the hardness of the 
process. 

3 Heavy quark production in the /cy-factorization 
approach 

The transverse momenta of the incident gluons in the small-x region result in the kr- 
factorization approach from a^lnfcy gluon diffusion. The diffusion is described by the 
so-called unintegrating gluon distribution - the function (^(x, q^) giving the gluon density 
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at the fixed fraction of the longitudinal momentum of the initial hadron, x, and the gluon 
virtuality, q^. These distributions should be found with the help of evolution equation 
and the experimental data. Unfortunately, it is not yet fulfilled. At very low x, that 
is to leading log(l/x) accuracy, the unintegrating gluon distribution are approximately 
determined via the derivative of the usual structure function: 

4v/2.'*^M!)1. (12) 

Such a definition of (^(x, q^) enables to treat correctly the effects arising from the gluon 
virtualities. 

Although generally is a function of three variables, x, and g^, the transverse 
momentum dependence is comparatively weak since ~ —q^ for small x in LLA in 
agreement with g^-dependences of structure function. Note that due to QCD scaling 
violation at larger g^, the value of ^p{x, q^) for the realistic structure functions increases 
faster with decreasing of x. Therefore the larger qx becomes important at small x in 
the numerical calculations. 

The exact expression of the gluon distribution can be obtained as a solution of 
the evolution equation which, contrary to the parton model case, is nonlinear due to 
interactions between the partons in the small x region. The calculations of the qx 



gluon distribution in leading order using the BFKL theory p4l result in differences from 



our (^(x, q ) function given by Eq. (12) by only about 10-15%. 



Here we deal with the matrix element accounting for the gluon virtualities and 
polarizations. Since it is much more complicated than in the parton model we consider 
only the LO of the subprocess gg — *■ QQ which gives the main contribution to the heavy 
quark production cross section at small x, see the diagrams a, b and c in Fig. 6. The 
lower and upper ladder blocks present the two-dimensional gluon functions ^{xi^ qf) and 
¥?(x2,g|). 

Strictly speaking, Eq. (12) may be justified in the leading log(l/x) limit only. The 
unintegrated parton distribution /a(x, qr, /i) determines the probability to find a parton 
a initiating the hard process with the transverse momentum qx (and with factorization 
scale /i). To restore the function /a(x, qx, /u) on the basis of the conventional (integrated) 
parton density a(x, A^) we have to consider the DGLAP evolution^ 



da a. 



(91nA2 27r 



Paa' {Z)a'{-, \')dz - a(x, A^) Y: Pa'a{z')dz' 



(13) 



(here a(x, A^) denotes xg{x, A^) or xg(x, A^) and Paa'{z) are the splitting functions). 

The first term on the right-hand side of Eq. (13) describes the number of partons 6a 
emitted in the interval A^ < gf^ < A^ + (5A^, while the second (virtual) term reflects the 
fact that the parton a disappears after the splitting. 

^For the g ^ gg splitting we need to insert a factor z' in the last integral of Eq. (13) to account for 
the identity of the produced gluons. 
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o b c d 



Fig. 6. Low order QCD diagrams for heavy quark production in pp (pp) collisions via 
gluon-gluon fusion (a-c) and the diagram (d) formally violating the factorization and 

being necessary for gaude invariance. 
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The second contribution may be resumed as giving survival probability Ta that the 
parton a with the transverse momentum qx remains untouched in the evolution up to 
the factorization scale 



Ta{qT,^) = exp 



■ , -i; 1^ / Pa'a{z )dz 



(14) 



For the case of one loop QCD running coupling as{pT) = 47r/61np|,/A^ the factor 
Ta{qT,fJ^) can be written down explicitely. In particular, the gluon survival probability 
(which enters our formulae) reads: 

TaidT, /.) = 1^:^^ ■ exp {M. [in ^ - In f In ^ - 2i?,(g^, ;.) + 3/2E,{qx, ;.)- 



-2/3E3(gr,/i) + l/4^4(gT,/i)] + 



(15) 



+^rzF [3Ei{qT, /i) - 3E2(gT, /i) + ^E^igr, 



where 



EkiqT,IJ') 



'A' 



[Et{k\n{p/A) - Ei{k\n{qT/A))] 



(16) 



and the integral exponent Ei(z) = /f ^ y expt; np and Nc are the number of light quark 
flavoures and the number of colours, respectively, and b = 11 — ^np. 



Thus the unintegrated distribution fa{x,qT,fi) reads 

»1-A 



fa{x,qT,ll) 



27r Jx 



Paa'{z)a' {^,q^] dz 



Ta{qT,fJ') 



(17) 



where the cut-off A = qr/fJ' is used both in Eqs. (14) and (17) |]45| , ^ 



We have to emphasize that fa{x,qT, fJ') is just the quantity which enters into the 
Feynman diagrams. The distributions fa{x,qT,fJ') involve two hard scale^: gr and the 
scale n of the probe. The scale n plays a dual role. On the one hand side it acts as 
the factorization scale, while on the other hand side it controls the angular ordering of 
the partons emitted in the evolution |4^, ^ The factorization scale /i separates 



the partons associated with the emission off both the beam and target protons (in pp 
collisions) and off the hard subprocess. For example, it separates emissions off the 
beam (with polar angle 6 < 90° in c.s.m.) from those off the target (with 6 > 90° in 
c.s.m.), and from the intermediate partons from the hard subprocess. This separation 
was proved in |^ ^ ^ and originates from the destructive interference of the different 
emission amplitudes (Feynman diagrams) in the angular boundary regions. 

If the longitudinal momentum fraction is fixed by the hard subprocess, then the limits 
of the angles can be expresseed in terms of the factorization scale /i which corresponds 
to the upper limit of the allowed value of the s-channel parton kp- 



^This property is hidden in the conventional parton distributions as qx is integrated up to the scale 
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Since the parton distributions contain two scales, one has to deal with a complicated 
evolution equation for fa{x,qT, fi)- On the other hand since the evolution is controlled 
by the emission angle only, the distribution /^(x, gr,/^) can be obtained from a single 
scale evolution equation, as it was done in Eq. (17). 

In the leading log(l/x) (i.e. BFKL) limit the virtual DGLAP contribution is ne- 
glected. So Tq = 1 and one comes back to Eq. (12) 

/r^"(^,gT,/i) = ^^^ ,x = gT. (18) 

In the double log limit Eq. (17) can be written in the form 

/i^^^(x,gT,/i) = \aix,X')TaiX,fi)]^_ , (19) 



X=qT 



which was firstly proposed by ||SD|. In this limit the derivative dTa/d In X"^ cancels the 



second term of the r.h.s. of Eq. (13) (see |^ for a more detailed discussion). 



Finally, the probability /a(x, gT,At) is related to the BFKL function ip{x,q^) by 

^{x,q^) = 4v^7rVa(x,gT,yu) ■ (20) 



Note that due to a virtual DGLAP contribution the derivative da{x,\^) / can 
be negative for not small enough x values. This shortcoming of Eq. (18) is partly 
overcome in the case of Eq. (19). Unfortunately the cut-off A used in a conventional 
DGLAP computation does not depend on the scale /x. To obtain an integrated parton 
distributions it is enough to put any small A ^ l[[ 

On the other hand, in the survival probability Eq. (14) we have to use the true 
(within the leading log approximation) value A = gT//^- Hence for a rather large qt 
(~ /i) and X even the DDT form Eq. (19) is not precise. Only the expression (17) with 
the same cut-off A in a real DGLAP contribution and in the survival probability (14) 
provides the positivity of the unintegrated probability /^(x, g^,/^) in the whole interval 
< X < 1. 

Of course, just by definition /a(x, gT,At) = when the transverse momentum qx 
becomes larger than the factorization scale /i. 

In what follows we use the Sudakov decomposition for the quark momenta pi^2 

through the momenta of colliding hadrons pa and Pb {p\ = P% — 0) and the trans- 
verse momenta pi,2T'- 

Pl,2 = Xi^2PB + yi,2PA + Pl,2T- (21) 

"^There is a cancellation between the real and virtual soft gluon DL contributions in the DGLAP 
equation, written for the integrated partons (including all fcr < m) • The emission of a soft gluon with 
momentum fraction (1 — z) — > does not affect the x-distribution of parent partons. Thus the virtual 
and real contributions originated from 1/(1 — z) singularity of the splitting function P{z) cancel each 
other. 

On the contrary, in the unintegrated case the emission of soft gluon (with > kx) alters the 
transverse momentum of parent (t-channel) parton. Eq. (17) includes this effect through the derivative 
^T{kj,,^i^)/^kl. 
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The differential cross section of heavy quark hadroproduction has the following form:^ 



'pp 



1 1 f 

'^Js)^ J ^^^ir^^^2T5(giT + q2T - Pit - P2t) 



dyldy*2(PpiTd?p2T (27r)8 (s) 

X ^^^{qly)^{qlx)\MQQ\\ (22) 

Here s = 2pxPB , Q'i,2T are the gluon transverse momenta and yl 2 are the quark rapidities 
in the hadron-hadron c.m.s. frame, 

xi= ^ e-y'^, X2= ^ e-y*^, x = Xi+X2 



y^=^ey^, y^ = ^ey^, y = yi + y2 (23) 

|MqqP is the square of the matrix element for the heavy quark pair hadroproduction. 



Contrary to the mention of p4|, the transformation Jacobian from x,y to yl,y2 is 



accounted in our matrix element. 
In LLA kinematic 



qi ^ ypA + qiT, q2 ^ + q2T- (24) 



so 

qf - -qlr^ ql - -qh- (25) 

2 2 

(The more accurate relations are q\ = ~"[~; 92 — t)ut we are working in the 

kinematics where x,y <^ 1). 

The matrix element M is calculated in the Born approximation of QCD without 
standard simplifications of the parton model. 

In the axial gauge p^A^ = the gluon propagator takes the form D^^^q) = d^^{q)/q'^, 
d,M) = - Wb + q''pi)/iPBq) . (26) 

For the gluons in t— channel the main contribution comes from the so called "non- 
sense" polarization g^^^, which can be picked out by decomposing the numerator into 
longitudinal and transverse parts: 

= 2(p>:; + p^M/s + ^ 2p>:;/s = g;^ . (27) 

The other contributions are suppressed by the powers of s. It is easy to check that in 
axial gauge (pB]) p^d^^ = and p^d^^ ~ —q^/y- Thus we effectively get 

d,M ^ -2^. (28) 
ys 



^We put the argument of 015 to be equal to gluon virtuality, which is very close to the BLM 
scheme (see also pl|). 
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Another way to obtain the same result is to use the transversahty condition. Since the 
sum of the diagram Fig. 6a-c is gauge invariant (see below) the product gi^iM^ = 0, here 
denotes the amplitude of the gluon qi and hadron pb interaction described by the 
lower part of Fig. 6a-c. Using the expansion p^ ) for the qi momentum we obtain 



y 

which leads to Eq. (pSf) or 

d,M) - 2^, (29) 
xys 

if we do such a trick for the vector pb too 0. Both these equations for d^i, can be used 



but for the form (28) one has to modify the gluon vertex slightly (to account for several 



ways of gluon emission - see ref. p^ ): 



^eff 



xys 



{xys - qlrr) qiT - qIt^At + 2x (girg2T) Pb ■ (30) 



As a result the colliding gluons can be treated as aligned ones and their polarization 
vectors are directed along the transverse momenta. Ultimately, the nontrivial azimuthal 
correlations must arise between the transverse momenta piT and p2T of the heavy quarks. 

From the formal point of view there is a danger to loose the gauge invariance in 
dealing with the off mass shell gluons. Say, in the covariant Feynman gauge the new 
graphs (similar to the " bremsstruhlung" from the initial or final quark line, as it is 
shown in Fig. 6d) may contribute in the central plato rapidity region. However this is 
not the fact. Within the "semihard" accuracy, when the function Lp{x,q^) collects the 
terms of the form ag(lng^)"(ln(l/x))™ with n + m > k, the triple gluon vertex ( pOD 
includes effectively all the leading logarithmic contributions of the Fig. 6d type ||4^, 
For instance, the upper part of the graph shown in Fig. 6d corresponds in terms of the 
BFKL equation to the t-channel gluon reggeization. Thus the final expression is gauge 
invariant (except a small, non- logarithmic, 0{as) corrections) 0. 

Although the situation considered here seems to be quite opposite to the parton 
model there is a certain limit |2^ , in which our formulae can be transformed into parton 
model ones, namely when the quark transverse momenta, pi,2r, are much larger than 
the gluon ones, gi,2T- 



^Having in mind this trick one can say that the matrix element is gauge-invariant in the kj'- 
factorization approach. However the polarization vectors of incoming gluons qi and q2 are not arbitrary 
taken but fixed as —q^j,/y and —ql^rp/x, respectively (see Q for more detail). 

^'^Taking the gluon polarization vector as —qix/y '^^-'^ completely eliminate the leading logarithm 
terms arising from Fig. 6d. 
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4 Qualitative difference between /ct- factorization ap- 
proach and parton model 



Eq. (p2D enables to calculate straightforwardly all distributions concerning one-particle 
or pair production. To illustrate the difference between our approach and the conven- 
tional parton model we present first of all the results of calculations of charm production 
{rric = 1.4 GeV ^) with high pt at three energies, y/s = 1 TeV, 10 TeV and 10^ TeV 



and with the same value of 

= ^ = 0.02 , (31) 
i.e. Pt = 10 GeV, 100 GeV and 10^ Gev for the above energies. 



We will use QCD scales jj^f = rriT = yml + p^ and = rric, i.e. fixed coupling 
with iV/ = 3 and A = 248 MeV |30|. 



However there exists a problem coming from the infrared region, because the func- 
tions ipi^x^ql) and (p{y,qf) are unknown at the small values of g^2- Moreover, for the 
realistic gluon structure functions the value (f{x, q^) is positive in the small-x region and 
negative in the large-x region. The boundary between two regions, where i~p{x, q^) =0, 
depends on the g^, therefore the relative contributions of these regions is determined by 
the characteristic scale of the reaction, say, by the transverse momentum pt- 

To avoid this additional problem, we present the numerical calculations both in the 
/cT-factorization approach and in the LO parton model, using the "toy" gluon distribu- 
tion 

xG{x,q^) = (l-xflniq^/Ql) , (32) 

for g2 > Ql and xG{x, q^) = for q^ < Ql with = 1 GeV\ 

After the calculations of charm production cross sections using, say, the VEGAS code 
|]52[] (that is the optimized Monte Carlo program for integrating the multidimentional 
expressions) and with Eqs. (12), (22), (32) one can see, that the values da/dxx at xt = 
0.02 are about 4-5 times larger in the fc-r-factorization approach than in the LO parton 
model, Eq. (1). Really this difference should not be considered as very large because, 
as it was discussed, an essential part of NLO and NNLO corrections is already included 
in the /cy-factorization, and it is known that the sum of the LO and NLO contributions 



in the parton model is about 2-3 times larger than the LO contribution ||55[| . 

To show the most important values of the variables 5i,2t in Eq- (|2^) , as well as the 
kinematical region producing the main difference with the conventional parton model, 
we plot by dots in Fig. 7 the results of the /cT^-factorization approach with the restrictions 
|5'i,2t| < q-max, as a function of qmax- The running coupling is fixed as 0:5 (m^) instead 
of <ys{(li2) in Eq. (0). The dashed lines in Fig. 7 are the conventional parton model 



predictions with QCD scales fip = y^nl + Pr-, and ixr = rric. One can see that the 
/cr-factorization predictions exceed the parton model results when qmax > Pt- 

Let us check now that the fcT-factorization predictions really coincide with the parton 
model ones when the quark momenta, pi,2T, are much larger than the gluon momenta. 
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Fig. 7. The values of da/dxT for charm hadroproductfon at xt = 0.02 in 
/cT-factorization approach as a function of upper hmits of integrals in Eq. (22) 
(points); the conventional parton model values (dashed lines) and the values of 
right-hand side of Eq. (34) (solid curves). 
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Fig. 8. The values of da / (Ixt/ dhiqmax for charm hadropro duct ion at xt = 0.02 in 
/cr-factorization approach as a function of upper hmits of integrals in Eq. (22) (points). 



?i,2T [iH- However we have to compare them at the same values of structure functions. 
When we submit Eq. ( ^2]) to the conditions |gi,2T| < <imax and neglect the gi,2r depen- 
dence in the matrix element, we recover the parton model approximation, and get the 
result proportional to xG{x, q^ax) 'yG{y, Qmax) due to the direct consequence of Eq. ( [l2l) 
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xG{x,q'^) = xG{x,Ql) + 



^{x,ql)dql 



(33) 



At the same time the original parton model yields the result proportional to xG{x, /i^^)- 
yG{y, fi'jp) with fip = \Jm1 + pj^. Hence we expect the parton model to coincide with 
our calculations for the gluon distribution Eq. (|32D , px ^ "^c and |gi,2T| < Imax after 
multiplying by an appropriate factor: 



da 

dxr 



da 



dxi 



PM 



Hp't/Q'o) 



(34) 



The right-hand side of Eq. (|3^ ) presented in Fig. 7 by sohd curves shows a good 
agreement with fc^-factorization approach (open dots) even when qmax is only slightly 
smaller (at the highest energy) than pt- The same values da/dxr, as in Fig. 7, but dif- 
ferentiated with respect to In qmax are presented in Fig. 8 for fcr-factorization approach. 
It exhibits, especially at the high energies, the logarithmic growth with qmax, until the 
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Fig. 9. The diagrams which are important in the case of one-particle distributions of 
heavy quark with large pt- High pt flows are shown by arrows. The situation similar 
to the LO parton model (a). The case, possible in NLO parton model, large pt of the 
quark is compensated by hard gluon and the fermion propagator is near to mass shell 

(b). The numerically small contribution, when large transverse momenta of heavy 
quarks are compensated by two hard gluons, whereas the fermion propagator is near to 

mass shell (c). 



value Qmax ~ Pt- There is the narrow peak in this region, which provides about 70-80 
% of the da/dxT cross section integrated over the whole gi_2T phase space. Its physical 
nature seems to be quite clear. There are two kinematical regions for t-channel and 
M-channel diagrams. Fig. 6 a,b, giving comparatively large contributions to da/dxT for 
the high energies and high pt heavy quark production. One of them comes from the con- 
ventional parton model kinematics when both initial gluon transverse momenta, qi^2T-, 
are very small compared to Pi^2Ti see Fig. 9a. Another large contribution appears from 
the region where, say, qix ~ Pit, whereas g2T and p2T are comparatively small, as it is 
shown in Fig. 9b. In this case the quark propagator, — — "^q) = 1/(92— P2 — '^q), 
is close to the mass shell and gives rise to the narrow peak shown in Fig. 8. The general 
smallness of such processes comes from high-virtuality gluon propagator in Fig. 9b, and 
it is of the same order as in the case of Fig. 9a, despite the diagram Fig. 9a is suppressed 
by the fermion propagator for the large rapidity difference between the two produced 
heavy quarks. 

The diagram shown in Fig. 9b can be considered ^ as one of the most important 
contribution to the NLO parton model in the high energy limit, because of spin-one 
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Table 1: The ratios of cc pair production in fc-r-factorization approach and in LO parton 
model 



yi, TeV 


0.3 


1 


10 


100 


1000 


Rtot 


4.0 


4.0 


4.0 


3.9 


3.9 


R{xt = 0.02) 


3.4 


4.5 


5.5 


5.4 


5.2 



gluon exchange in the t-channel. Due to these factors, combinatorics and interference 
between diagrams, the resuhing /cy-factorization contribution to the peak at qt ~ Pt 
in the total da/dxT value in Fig. 8 is several times larger than the LO parton model 
contribution. In the Table 1 we present the calculated ratios of the total heavy quark pair 
production cross section, Rtot, and the ratio of da/dxT, R{xt) at xt = 0.02 (integrated 
over rapidity). 

One can see that at fixed xt the relative contribution of the discussed peaks at first 
increases due to increase of the phase space. This contribution is saturated at the high 
enough energy {^/s ~ 10 TeV). After this the relative contribution of the LO parton 
model increases logarithmically with pt, and have to dominate at the extremely high 
energies and pt, in academical asymptotic. 

These results are confirmed by those presented in Fig. 10, where we reproduce the 
data from Fig. 7 for LO parton model and /cT-factorization approach with the condition 
|?i,2t| < Qmax- FoT comparisou we show by stars the /cj-factorization predictions for all 
values qiT with the only condition |g2T| < Qmax, versus qmax- The last results are above 
the LO parton model even at not too large q2T because the peaks, discussed above, are 
included into the integral over qiT- 

Let us note that the calculation of da/dxT at ^/s = 10 TeV, xt = 0.02, and with 
restriction |gi,2T| > Pt/2 (see Fig. 9c) gives only about 2% of the total value of da/dxx- 

The essential values of gi,2T increase in our calculations when the transverse momen- 
tum, Pt of the detected c-quark increases. At the high initial energy gi 2t ~ Pt- In the 
kick language it means that the {k^) value also increases. 



5 Total cross sections and one-particle distributions 

Eq. (22) enables us to calculate straightforwardly all distributions concerning one- 
particle or pair production. One-particle calculations as well as correlations between 



two produced heavy quarks can be easily done using, say, the VEGAS code . 

However there exists a principle problem coming from the infrared region. Since 
the functions ip{x, g|) and (/?(?/, q\) are unknown at small values of g| and , i.e. in 
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Fig. 10. The same as in Fig. 7 (points and dashed curves) together with the values of 
da/dxT (stars) for the case when only g2T upper limit integration is bounded by Qmax- 



nonperturbative domain we use Eq. (33) and rewrite the integrals in the Eq. (22) as 

[ d^qiTd'^q2T5{qiT + q2T - Pit - P2T)^^^^^^^fg{y:qiT: fl)fg{x,q2T: ^)\MQQf = 

J qi ?2 

= {4V2n'a,{ml)f xG{x, Ql) yG{y, Ql) T\QI f^') 0-^^] + (35) 

V 1112 / 912^0 

+ 4V2TT^as{ml)xG{x, Ql) T{Ql, /x^) / dq^j- S{qiT - Pit - P2t) x 

■'Ql 

91 \ 92 



92^0 



CXD 



+ 4^2 7r2a,(m^) yG{y, Ql) T{Ql, /x") / dqi^ 5{q2T - Pit - P2t) x 



X % fa{x,q2T,IJ') 



/•CXD /"OO 

+ / d'^qiT / , d^q2T5{qiT + 92T - Pit - P2t) x 

X '^^h{y,q.T,^^)h{x,q2T,^^)\MQQ\' , 
9i 92 
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where the unintegrated gluon distributions fg{x,qT,fi) are taken from Eqs. (14) and 
(17). In the numerical calculations we use the values /i^ = s and /x^ = s/4. 

The first contribution in Eq. (35) with the matrix element averaged over directions 
of the two-dimensional vectors Qit and g2T is exactly the same as the conventional LO 
parton model expression. It is multiplied by the 'survival' probability T^((5o,/W^) not to 
have transverse momenta qiT, q2t > Qo- The QCD scales are = mf. and = Ql 
(the same value as in Eq. (33), we assume Ql = 1 GeV^) The sum of the produced 
heavy quark momenta is taken to be exactly zero here. 

The next three terms contain the corrections to the parton model due to gluon po- 
larizations, virtualities and transverse momenta in the matrix element. The relative 
contribution of the corrections strongly depends on the initial energy. If it is not high 
enough, the first term in Eq. (35) dominates, and all results are similar to the con- 
ventional LO parton model predictions. In the case of very high energy the opposite 
situation takes place, the first term in Eq. (35) can be considered as a small correction 
and our results differ from the conventional ones. It is necessary to note that the abso- 
lute as well as the relative magnitudes of all the pieces in Eq. (35) strongly depend on 
the T-factor (i.e., when we use Eq. (17) or Eq. (12)). 

Before the numerical comparison it is necessary to remind that the NLO parton 
model actually results only in a normalization factor in the case of one-particle distri- 
butions, the shapes of LO and LO-I-NLO distributions are almost the same |^, |5|, |5l 



see Fig. 2. This means that we can calculate the K-factor Eq. (7), say, from the results 
for the total production cross sections, and restrict ourselves only to the LO calculations 
of pt or rapidity distributions multiplying them by the K-factors. 

The numerical values of the K-factors depend ||55| on the structure functions used. 



quark masses, QCD scales and the initial energy, the renormalization scale fiR being 
especially important. This seems to be evident, because the LO contribution is propor- 
tional to a|, whereas the NLO contribution is proportional to a|. However the structure 
of Eq. (2) is more important at the high energies, when small p values dominate. At 
p ^ the functions /^^ and /^^^ have the constant limits fg]^{p — > 0) ~ 0.1 and 
fS^P ^ 0) ~ -0.04, while ffj{p ^ 0) ^ 0. Therefore due to Eq. (2) the K factor 
value at the high energies is mainly determined by the ratio p/rriQ. 

First of all let us present the role of the T-factors, Eq. (14). In Fig. 11 we show 
the ratios of the last term of Eq. (35) with the T-factors in both gluon distributions 
fg{x, qiT, p) to the same values calculated for Tg{x, qiT, p) = 1 in one unintegrated gluon 
distribution. The ratios are plotted as functions of qiT for the cases of charm and beauty 
production at ^/s = 14 TeV and /x^ = s. The heavy quark transverse momenta are fixed 
to be 20 GeV/c. In both cases the factors Ta{qu) are rather small at small qiT and 
Tg{x, qiT, p) 1 at > gir- 

Now let us compare the numerical results predicted by the parton model and by the 
fcr-factorization approach. 

The energy dependences of the total cross sections of cc and bb pair production 
are presented in Fig. 12. As was mentioned, at comparatively small energies the first 
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Fig. 11. The role of T-factor, Eq. (14) in the calculation of charm (solid curve) and 
beauty (dashed curve) production with px = 20 GeV/c at ^/s = 14 TeV and /i^ = s. 



term in Eq. (35) dominates and the results of the /c-r-factorization approach should be 
close to the LO parton model prediction. Actually the first results are even smaller 
due to the presence of the T-factor in Eq. (17). However the /cT-factorization approach 
yeilds a stronger energy dependence than the LO parton model both for cc and bb 
production. This can be explained by additional contributions appearing at very high 
energies in the /cr-factorization approach when the transverse momentum of a heavy 



quark is compensated by the nearest gluon, see |^| 

One-particle pr distributions, da/dpT, calculated in the /c-r-factorization approach 
and in the LO parton model are presented in Fig. 13. In all cases the fc^-factorization 
predicts broader distributions. The average values of pt of the produced heavy quarks 
are rather different in these two approaches, as one can see from Table 2. 

Table 2 The average values of charm and beauty quark transverse momenta <Pt> 
(in GeV/c) in the fc^-factorization approach with fi'^ = s and in the LO parton model. 
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Fig. 12. The total cross sections for charm and beauty hadroproduction in the 
/cT-factorization approach with unintegrated gluon distribution fg{x,qT,fJ^) given by 
Eq. (17), for //^ values in Eq. (14) equal to s (sohd curves), and s/4 (dash-dotted 
curves) , and in the LO parton model (dashed curves) . 
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Fig. 13. pT-distributions of c-quarks (a) and 6-quarks (b) produced at different 
energies. Dashed curves are the results of the LO parton model. Solid curves are 
calculated with the unintegrated gluon distribution /^(a;, g^,/^) given by Eq. (17), for 
/i^ values in Eq. (14) equal to s and dash-dotted curves are calculated for /i^ = s/4. 



This seems to be very natural, because, contrary to the LO parton model, a large 
Pt of the one heavy quark can be compensated not only by the pt of the other heavy 
quark but also by neighbour gluons emission. Fig. 9b. 

The rapidity distributions of produced heavy quarks presented in Fig. 14 show that 
the main part of the difference between the /cj^-factorization approach and the LO parton 
model comes from the central region. 



6 Two-particle correlations 

We saw from the previous section that there is not so large difference in our results for the 
total cross sections and one-particle distributions obtained in the /cT-factorization and 
in the LO parton model. The predictions of the NLO parton model for these quantities 
differ from the LO parton model only by a normalization factor of 2-2.5 [§, ^, ^, |54 |. 



Hence the difference between our predictions and the NLO parton model should be 
small. 

The calculations of two-particle correlations in different approaches are more infor- 
mative. The simplest and/or most informative correlation here is the distribution of the 
total transverse momentum of the produced heavy quark pair, Ppair- In the LO parton 
model Ppair = Pit + P2T = Qit + g2T, therefore if qiT = q2T = 0, then da/dppair is a 
^-function with Ppair = 0. Thus the Ppair distributions provide the direct information 
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Fig. 14. Rapidity distributions of c-quarks (a) and 6-quarks (b) produced at different 

energies. Dashed curves are the results of the LO parton model. Solid curves are 
calculated with the unintegrated gluon distribution fg{x,qT,fi) given by Eqs. (17), for 
fj? values in Eq. (14) equal to s and dash-dotted curves are calculated for /i^ = s/4. 



about the transverse momentum distribution of the incident partons. 

It is clear that if Qit <C Pit, then the distributions in Ppair should be narrower in 
comparison with the one-particle pt distributions. In this case the Weizsaecker- Williams 
approximation should be valid and one can believe that the parton model reflects the 
real dynamics of the interaction. In the opposite case, qix ~ PiT, the large transverse 
momentum of the produced heavy quark can be compensated not by the other quark, 
but by a high-p^ gluon. We have shown (see Sect. 4), that about 70-80% of the 



total cross section of high-p^ quark production at high energies originates from such 
processes, when the heavy quark propagator is close to the mass shell. 

We calculate the values of da/dppair for charm (a) and beauty (b) production in the 
/cT-factorization approach using the unintegrated gluon distribution Eqs. (14), (17) with 
scale values /x^ = s and /i^ = s/4 (the last value only for -y/i = 14 Tev). Our results for 
pair production at different initial energies are shown by solid curves in Fig. 15. For the 
comparison we present by dashed curves the one-particle pr-distributions taken from 
Fig. 13, obtained in the same fc-r-factorization approach and with the same T-factor. As 
we put Ql = 1 GeV^ in Eq. (35), we can not distinguish between the initial gluons with 
qT equal to, say, 0.1 GeV/c and 0.9 GeV/c, so our first bin in the da/dppair distribution 
has the width 2 GeV/c which explains some irregular behaviour of the solid curves at 
the small pt- Naturally, all the solid and dashed curves are equally normalized at the 
same energy. 

At comparatively small energies, ^/s = 39 GeV and even at ^/s =630 GeV the 
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Fig. 15. The distributions of the total transverse momentum "Ppair for c-quarks (a) and 

6-quarks (b) produced at different energies (solid curves), calculated with the 
unintegrated gluon distribution fg{x,qT,IJ') given by Eq. (14) and (17), for /i^ values in 
Eq. (14) equal to s. Dashed curves show the one-particle (single) pr-distributions with 
the same /i^ taken from Fig. 13. Dash-dotted and dotted curves are the same 
calculations for pair and single production at y/s = 14 TeV with /^^ = s/4. 
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Fig. 16. The calculated distributions of the rapidity difference between two c-quarks 
(a) or 6-quarks (b) produced at different energies in the kr factorisation approach, 
calculated with unintegrated gluon distribution fg{x,qT,lJ') given by Eq. (17), for //^ 
values in Eq. (14) equal to s (solid curves) and ji"^ — s/A (dash-dotted curves). Dashed 

curves show the LO PM predictioms. 



distributions da/dppair are narrower than the one-particle distributions da/dpr- This 
means that the transverse momenta of the produced heavy quarks almost completely 
compensate each other. However the situation changes drastically with increasing of 
the initial energy. Starting from comparatively small pt, the difference between the 
curves decreases with energy. At y/s = 14 TeV the distributions are similar both in 
the cases of cc and bb production. This means that the production mechanism changes 
in the discussed energy region. At ^/s = 14 TeV the transverse momentum of the 
produced heavy quark is balanced more probably by one or several gluons, because 
the contribution with large virtuality in the quark propagator is more suppressed in 
comparison with the large virtuality in the gluon propagator. 

The discussed behaviour depends on the value of the scale /i^ in the T-factor, Eq. (14). 
The similar calculation at energy ^/s = 14 TeV with fi^ = s/4 is shown in Fig. 15 for 
pair and single production by dash-dotted and dotted curves, respectively. Here the 
difference between these two curves is more significant and becomes larger for lower 
energies. 

The distributions of the produced heavy quark pair as a function of the rapidity gap 
— Iuq ~ VqI between quarks are presented in Figs. 16. Difference between the LO 

PM and the /cr-factorization predictions is not large again except for the region of very 

large Ay. 

Another interesting correlation is the distribution of the azimuthal angle (j). Pre- 
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liminary results for the azimuthal correlations in the fcr-factorization approach were 



considered in p3[. The main difference in the information coming from da/dppair and 
da/d(j) distributions is due to the comparatively slow heavy quark. It gives a negligi- 
bly small contribution to the da/dppair in the first case, whereas in the second one each 
quark contributes to the distribution da/d(j) practically independently of its momentum, 
so all corrections coming from quark confinement, hadronization and resonance decay 
can be important. 

As was discussed above, the first contribution in Eq. (35) is the same as the con- 
ventional LO parton model in which the angle between the produced heavy quarks is 
always 180°. However the angle between two heavy hadrons can be slightly different 
from this value due to hadronization processes. To take this into account we assume 
that in this first contribution the probability to find a hadron pair with azimuthal angle 
(j) = 180° — 01 is determined by the expression 

wi{(pi) = I f ^ , (36) 
^ PI + Pt 

where ph = 0.2 Gev/c is a transverse momentum in the azimuthal plane coming from 
the hadronization process. The other contributions of Eq. (35) result in a more or less 
broad distribution so we neglect their small modification due to hadronization. 

The /cT-factorization approach predictions for the azimuthal correlation of heavy 
quarks produced in pp collisions are presented in Fig. 17. One can see that they change 
drastically when the initial energy increases from fixed target to the collider region. 

There exists a lot of various correlations which have never been considered both the- 
oretically and experimentally. Let us discuss several of them. In Fig. 18 we present the 
distribution on heavy quark pair transverse momentum in events, when the transverse 
momentum of one quark is fixed in the interval piT = 19-21 GeV/ c. These distributions 
are very sensitive to the scale value fj? in Eqs. (14), (17). In the case of /i^ = s both 
the charm and beauty distributions are practically flat for ppair < Pit (with evident 
exception for the kinematical minimum at very small Ppair) and decrease rather fast for 
Ppair > Pit- Such a behaviour again shows that the high transverse momentum of one 
heavy quark can be compensated by the second quark (region of comparatively small 
Ppair as well as by the hard gluon (region of Ppair ~ Pit)- For the smaller scale s/4 in 
Eqs. (14), (17) the emission of hard gluons is suppressed and the dashed curves decrease 
immediately after the kinematical minima. Thus we can conclude that the measurement 
of this distribution allows to find the most reasonable scale value. 

The scale problem can be solved if we will consider some distributions in different 
azimuthal angle regions. For example, the one-dimentional transverse momentum dis- 
tribution is insensitive to the scale for back-to-back production (A0 = 120° — 180°) and 
rather sensitive for the smaller azimuthal angles, see Figs. 19 and 20. 

The distribution over Ay at LHC energy is predicted to have minimum at the small 
azimuthal angles and Ay values in the case of comparatively small scale //^ = s/4 and 
to be more fiat for the larger scale, see Figs. 21, 22. 
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Fig. 17. The calculated azimuthal correlations of charm (a) and beauty (b) pair 
production for fi"^ = s (solid curves) and fi'^ = s/i (dash-dotted curves) at the energies 
equal to y/s = 14 TeV (upper curves at small 0), 1.8 TeV, 630 GeV and 39 GeV (the 

latter one only for charm production). 



7 Conclusion 

We have compared the conventional LO Parton Model (PM) and the fc-r-factorization 
approach for heavy quark hadroproduction at fixed target and collider energies using 
both "toy" and realistic gluon distributions. Both the transverse momenta and rapidity 
distributions have been considered, as well as several types of two-particle correlations, 
such as the distribution of the rapidity gap between two heavy quarks, their azimuthal 
correlations and distributions of the total transverse momentum of the produced heavy 
quark pair, Ppair- 



It has been shown in Sect. 4 that the contribution of the domain with strong 



qt ordering (gir,Q'2T = y^Q + Pt) coincides in the /ct- factorization approach 

with the LO PM prediction. Besides this a numerically large contribution appears at 
high energies in /cT-factorization approach in the region gi,2r > ^t- This kinematically 
relates to the events where the transverse momentum of heavy quark Q is balanced not 
by the momentum of antiquark Q but by the momentum of the nearest gluon. 

This configuration is associated with the NLO (or even NNLO, if both qix, q2T > tt^t) 
corrections in terms of the PM with fixed number of flavours, i.e. without the heavy 
quarks in the evolution. Indeed, as was mentioned in [|ll, up to 80% of the whole NLO 
cross section originates from the events where the heavy quark transverse momentum 
is balanced by the nearest gluon jet. Thus the large "NLO" contribution, especially 
at large is explained by the fact that the virtuality of the t-channel (or -u-channel) 
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Fig. 18. The distributions of the total transverse momentum Ppair for c-quarks (a) and 
6-quarks (b) produced at different energies (sohd curves) and with the transverse 
momentum of a one heavy quark restricted to the interval 19-21 GeV/c, calculated 
with the unintegrated gluon distribution fg{x,qT,IJ>) given by Eq. (14) and (17), for /x^ 
values in Eq. (14) equal to s. Dashed curves show the same distributions with scale 

= s/4. 




Fig. 19. pT-distributions of c-quarks in different azimuthal angles produced at different 

energies. 
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Fig. 20. pT-distributions of 6-quarks in different azimuthal angles produced at different 

energies. 




Fig. 21. Tfie calculated distributions of the rapidity difference between two c-quarks 
produced at different azimuthal angles in the kr factorisation approach. 
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Fig. 22. The calculated distributions of the rapidity difference between two 6-quarks 
produced at different azimuthal angles in the fcy factorisation approach. 



quark becomes small in the region qt — Pt, and the singularity of the quark propagator 
^/{p ~ Q ~ "^q) in the "hard" QCD matrix element, M{qiT, ?2T,Pit,P2t), reveals itself. 

The double logarithmic Sudakov-type form factor T in the definition of the unin- 
tegrated parton density Eq (17) comprises an important part of the virtual loop NLO 
(with respect to the PM) corrections. Thus we demonstrate that the /^T-factorization 
approach collects already at the LO the major part of the contributions which play the 
role of the NLO (and even NNLO) corrections to the conventional PM. Therefore we 
hope that a higher order (in as) correction to the /cT-factorization could be rather small. 

Another advantage of this approach is that a non-zero transverse momentum of QQ- 
system {pTpair = Pit + P2T — Qit + 92t) is naturally achieved in the /c^-factorization. 
We have calculated the pTpair distribution and compared it with the single quark pt 
spectrum. At the low energies the typical values of pxpau- are much lesser than the 
heavy quark pt in accordance with coUinear approximation. However for LHC energy 
both spectra become close to each other indicating that the transverse momentum of 
second heavy quark is relatively small. The typical value of this momentum {pTpair = kr- 
kick) depends on the parton structure functions/densities. It increases with the initial 
energy (/c^^-kick increases with the decreasing of the momentum fractions x, y carried by 
the incoming partons) and with the transverse momenta of the heavy quarks, pt- Thus 
one gets a possibility to describe a non-trivial azimuthal correlation without introducing 
a large " phenomenological" intrinsic transverse momentum of the partons. 

It is necessary to note that the typical parton transverse momenta Qit and q2T 
increase in our calculations with the growth of the detected 6-quark momentum p^J*". 
In the language of /cy-kick it means that the values of {k^) also increase. 
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At the time being we would not like to compare our calculations with experimental 
dataQ. It is shown that in all cases when a realistic gluon distribution is used the 
difference between our approach and the conventional parton model is not large, so our 
results should be in agreement with the data up to, say, a factor 2 or 3. At the same 
time the disagreement of such an order seems to be senseless to discuss, because the 
unintegrated gluon distributions are not known with the needed accuracy (from the 
evolution equation and the experimental data). Instead of them we have used more or 
less reasonable approximation Eq. (17). Although it is qualitatively good it can lead to 
a numerical disagreement. In particular, the angular ordering |^ ^ implies that 
the cut-off A = Qt/ H in Eqs. (14), (17) should be replaced by A = qr/ifJ' + It)- From 
the formal point of view the difference is beyond the DGLAP LO accuracy. However 
with the new A = (3't/(/^ + (It) one gets the non-zero values of fa{x,qT,[Ji) even at the 
large qt > /i. This is especially important at the high energies where the essential 
values of x and z (in Eq. (17)) are very small. The contribution coming from qx > fi 
region could enhance the flux of colliding gluons up to a few times, and by this way 
explain the FNAL-Tevatron puzzle. These new data on the cross section of bb (or high- 
Pt prompt photon) production are 2-3 times larger than the conventional NLO PM 
QCD predictions |10|, At the moment we have no realistic unintegrated parton 
distributions which fit the data with accounting for the contribution from qx > fi- 

That is why the presented here results are treated mainly as the qualitatively ones, 
having the numerical accuracy on the level of factor 2-3. 

Another important point is that almost all presented results concern the heavy quarks 
rather than the hadron production. Of course, the hadronization leads to several im- 
portant effects, such as different yields (asymmetry) of D- and Z)-mesons production 
in pp collisions, see qualitative discussion in |5^, however their quantitative description 



needs additional phenomenological assumptions, see e.g. jS^, |5^, |5^, Q description 
of PYTHIA Monte Carlo code. 
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